Tailoring photon emission patterns in nanostructures 
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We investigate the photon emission in coupled quantum dots based on symmetry considerations. 
With the help of a new theorem we proved, we reveal the origin of the various emission patterns, 
which is the combinative symmetry in the time domain and spectrum domain. We are able to 
tailor the emission patterns and obtain emission spectra with odd harmonics only, even harmonics 
only, both odd and even harmonic components, or even the quenching of all components. These 
interesting emission patterns can be obtained in experiments by careful design of the nanostructures, 
which are of many applications in optical-electric nanodevices. 
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INTRODUCTION 

Photon emission in nanostructures plays a crucial role 
in modern electronic/optical devices. Generating high- 
order harmonics is one efficient up-conversion method 
for obtaining desired spectra (for example, THz spec- 
tra) from sources with lower frequencies. During the last 
several decades, much effort was made in the study of 
generation of high-order harmonics in atomic/molecular 
systems and nanostructures [ll-Q, as well as their appli- 
cations. For instance, Ahn et al recently proposed an ap- 
proach for THz wave generation by high-order harmonic 
wave based on semiconductor nanostructures driven by 
acoustic wave A scheme of electrically pumped 

photonic-crystal THz laser [H was developed by Chas- 
sagneux et al. A method for THz wave generation by 
Gigahertz wave was also suggested in 

Different harmonic components have interesting appli- 
cations. For example, even harmonics were used as a 
test wave to diagnose the fast time evolution of the cur- 
rent density Yet odd harmonics were often observed 
in experiments in atomic and molecular systems. The 
appearance of odd harmonics was attributed to the par- 
ticular inversion symmetry in central potential. There 
have been several studies on the dynamical symmetry Q 
and related selection rules @ in high-order harmonics 
generation. People have tried to get even harmonics by 
avoiding the selection rule [l0l - [l3 | . Some studies focused 
on the symmetry breaking for obtaining the even har- 
monics in molecules and atoms p^ - [l^ . The symmetry 
breaking may be realized in the molecule with two nu- 
cleus of different masses [l^. The even harmonics were 
found to appear in driven double quantum wells, where 
the potential was not of inversion symmetry. Also the 
radiation may occur at non-integer multiples of the fun- 
damental frequency. 

There have been many theoretical and experimental 
studies on high-order harmonics, and most of them fo- 
cused on the atomic and molecular systems [l5l - l27| . In 
spite of many studies on even harmonics generation, it 
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appears that the deep origin of different emission pat- 
terns is still unclear, and an effective way of generating 
emission spectra of specific pattern is lack. In this paper, 
we study the emission spectra of coupled quantum dots 
(QDs), and find the symmetry origin of various emis- 
sion patterns, as well as methods for generating different 
emission patterns including those with even harmonics. 
The main advantage of QD (artificial atom) and coupled 
quantum dots (CQDs) (artificial molecules) is their tun- 
ability. By carefully designing the material, the growth 
process, or adding appropriate gate voltages etc., one is 
able to adjust the energy levels/energy gaps in CQDs. 
One can further design the structure of CQDs (the rel- 
ative position of QDs, the inter-dot distances, the hight 
and width of the tunneling barriers) to tune the optical 
dipole between QDs. The optical coupling between QDs 
can also be changed by tuning the polarization of the in- 
cident light. Compared with one quantum dot with mul- 
tiple energy levels, the CQD system has the advantage of 
more tunability. By making full use of the tunability of 
coupled nanostructures, we propose an effective approach 
to tailor emission patterns. We reveal the origin of ap- 
pearance of odd/enen harmonics based on the symmetry 
considerations. It turns out that the emission pattern 
is determined not just by the inverse symmetry, but by 
a new type of symmetry, which is the combinative sym- 
metry in the time domain and spectrum domain. Based 
on our findings, we are able to obtain emission patterns 
with even harmonics only, odd harmonics only, or both 
even and odd harmonics, and even disappearance of all 
harmonic components. Our methods of generation of var- 
ious emission patterns in nanostructures have important 
applications. 



THEORETICAL FORMULISM 

We consider a CQD system with one energy level for 
each dot. The energy is Ei for the state |z) in the dot 
i, i = I...N. This CQD is driven by an external field 
E ~ F coa{u!ot)u, u unit vector. 

Under the dipole approximation, our system is de- 
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scribed by the Hamiltonian 

H = J2E^\^){^\+J2 cos(c.ot)N)01, (1) 

where Gij = Fu- flij are Rabi frequencies, flij = {i\er\j) 
the dipole between dot i and dot j. The equation of 
motion for the density matrix is written in the form (28| 



where the last term describes possible dissipation ef- 
fects (for instance that from the spontaneous phonon 
emission). We set ?i = 1 in the following. The time- 
dependent mean dipole moment can be calculated as 
P(t) = (er) = f^ijPijW- With the help of Fourier 
transformation we can obtain the emission spectrum 
S{u!) = I / dt cxp{iujt)P{t)\'^ . Now we present a theorem 
on the origin of the various emission patterns. Several 
examples will be given. 

Theorem For a quantum system described by the 
Hamiltonian (1), if there exists one symmetric operation 
Q, which is the time shift : t ^ t + T/2 (or -t + T/2, 
T = 27r/wo) combined with another operation in spa- 
tial/spectrum domain, i.e. Q — ■ 9 , such that the 
initial condition and the Hamiltonian are invariant (or 
H —H), and the dipole operator P has a definite 
parity, then the emission spectrum contains no odd/even 
component if operator P is even/odd. 

The proof goes as following: 

Proof We first consider the case: under the operation 
Q, H ^ H. The schrodinger equation = Hlip) 

remains invariant under Q transformation. If the initial 
condition remains unchanged (up to a phase), that is, the 
initial condition |V'(i = 0)) satisfies \'ip{t = 0)) = Q\tp{t = 
0)) = e-*T|^/'(i = 0)), then we haveJV'(i))^e-^T|V'(i))- 
Therefore P(i) = (i;{t)\P\^{t)) = {ij}{t)\P\^p{t)) = {^p{t + 
T/2)\n-^Pn\ij{t + T/2)) = {iP{t + T/2)\ ± P|V'(i + 
T/2)) = ±P{t + T/2), where for even P and "-" 
for odd P. Thus we have P(nwo) = / dte'^^^o^Pit) = 
±/dte"'^"(*-^/2)p(t) ^ ±(-l)"P(na;o). 

Similarly for the case: under the operation Q, H ^ 

-H. p{t) = {m\pm)) = {m\p\m) = m-t + 

T/2)\n-^Pn\tl;{-t + T/2)) = {^{-t + T/2)\±P\tlj{-t + 

T/2)) = im - T/2)\)* ± Pmt T/2)))* = m - 

T/2)|±P|V'(t-P/2)) = ±P{t-T/2), where we have used 
the fact that P{t) is real. Thus we also have P{nujo) ~ 

±(-l)"P(nL.o)- 

So we reach the results, 

P{nu!o) = for odd n, if P is even (under il); 

P{nu!o) = for even n, if P is odd (under 57). 

From above proof, one sees that the required ini- 
tial condition is \^p{t = 0)) = e^TQ|?/)(t = 0)) = 
e'^n\i;{T/2)) = e''^f7C/(r/2)|?/>(0)), where U is the time 
evolution operator. So the initial state needs to be an 
eigenvector of the operator VlU{T/2). This initial condi- 
tion is not very convenient for practical use since the 



time evolution operator is involved. Here we demon- 
strate that the initial condition could be replaced with 
\il}{t = 0)) = e*T2|-i/)(t — 0)), and the emission spectra 
have little change, as long as \Ei — Ej\/u]Q is small. 

Physically we can understand that for linear time evo- 
lution systems, very small change of initial condition does 
not lead to big change of the time evolution of the sys- 
tem. More precisely, from Floquet theorem, we have 
hAl^)) ~ Cae~''°*|wQ(t)), with \ua{t)) the time pe- 
riodic Floquet state, Cq, determined by the initial condi- 
tion. It is easy to sec that P(w) — ^, c*,Cq,A'q^q./ (a;), 

where Ka,a'{i^) = /dte*'^*e-'(^°-^°)*<,(i)Pua(<)- So 
for two different initial conditions with Ca,Za, such that 
Cq, — CqI < rj, r] a, sufficient small parameter, it is 

easy to show that |P(w) — ^'('^)| < J2a a' '^V^a.a', here 
P{uj), P{lu) are the emission spectra corresponding to the 
initial conditions with Cq, respectively. 

Now we show that the difference between the initial 
conditions satisfying |-0(t ~ 0)) ~ e'^Tl\ip{t — 0)) and 
V'(i = 0)) = e^^Q\^p{t = 0)) can be small enough, if 
E'i — E'j I/wq <C 1. The time evolution operator satisfies 
the equation idU/dt = HU with initial condition U{t = 
0) = I, I the identical operator. H = E + G cos{ujot)fi, 
where E = diag(Ei, En). The symmetric matrix 
ft can be diagonalized as A~^fiA ~ diag{\i, ...\n)- 
Using the transformation U = ABUA~^, B = 
diag(e-'^iG^W,...,e-*^«<^/(*)), f{t) = /„* dr cos(wot), 
one sees that idU/dt = EU, where E = B-^A~'^EAB, 
U{t = 0) = /. Then _it is clear that |J7(r/2) - I\ 
and \U{T/2) - /| = \A[U{T/2) - I]A-^\ can be smaU 
enough, if \E^ - Ej\/ujo is small. Since Q = W{T/2) = 
n + nU{T/2) -n = n + n{U{T/2) - /), standard pertur- 
bative calculation shows that the correction of initial con- 
dition due to perturbation Vl{U{T /2)—I) is small and one 
can use the initial condition \ilj{t = 0)) = e*Ti|V'(t = 0)) 
for the cases with \Ei — Ej\/ujQ In the above proof 

the dissipation hasn't been included, yet we would like to 
point out that if the dissipation is small, the main feature 
of emission spectra remains unchanged as verified by our 
numerical calculations shown in the next section. 

We have obtained one selection rule in the perturba- 
tion regime ( \Ei — Ej\/uj[) <C 1) and it depends on the 
initial condition. What may happen in the nonpertur- 
bation regime? In the regime with large \Ei — Ej\/uiQ 
and G/luq, most quasienergies are large and the corre- 
sponding quasienergy states in the expansion \ip{t)) = 
X]q Cae^^^"*-\ua{t)) lead to fast osciUation behavior, only 
the quasienergy state with the smallest quasienergy 
(modulo Tiujq) dominates, i.e., |V'(i)) = Cage~^'^''o^Ua„{t). 
If there exists one symmetric operation Q such that 
QHQ-^ = H, Qi-^Q~^ = i-^, then the state Q|uao) 
is also a quasienergy state with quasienergy If 
the quasienergy states are nondegenerate, Q|uqo) is just 
\uao) (up to a phase). Straightforward calculation shows 
that P{t) = -P{t + T/2) (Note that P is odd in this 
case). Therefore there is no even harmonics. One should 
notice that this selection rule is insensitive to the initial 
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condition. 

To have a clearer understanding of the physical pic- 
ture, we study the case of double dots (with two levels), 
since some explicit solution could be obtained. The prob- 
ability amplitudes /3,a for one electron in states |1) and 
2) satisfy the following equation 



Gcos{ujot) 
A 





where A — E2 ~ Ei. Define A = e" 
B = e^^'^^^^\a — (3). Then we have the equations for A 
and B 



(4) 



The time-dependent dipole P can be written as 

P{t) = 2ni2Re{a*l3} = ^ii2Re{A*A - B*B}/2. (5) 

After some algcbric calculation, we have the equation for 
P{t) 



^ = AIm{{a* + b*){a - b)e'^^f(''>} 
-A2/p*drcos[2G/(t) - 2G/(t)]P(t), 



(6) 



where a, b are the probability amplitudes for the elec- 
tron in the two dots at initial time t ^ 0. We have 
made no approximation and the above equation is exact. 
Here we see that the first term, i.e. the initial condi- 
tion is important for A <C 1. The dependence of emis- 
sion spectra on the initial states has also been noticed in 
j29| . P{t) becomes less sensitive to the initial condition 
for large A as we've discussed in the nonpcrturbation 
regime. For a particular type of initial condition a = b, 
i.e. pii{t = 0) = p22{t = 0) = 0.5, the first term on the 
right side of equation (6) vanishes. Then the equation 
can be solved for small A and one finds that there is no 
odd harmonic in the emission spectrum [3^. It agrees 
with our theorem as also shown in Fig. (lb). 



NUMERICAL RESULTS AND DISCUSSIONS 
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FIG. 1: Emission spectra in double dots with different initial 
conditions. A = £1 - £2 = 0.02, G12 = 17.0, and V = 
6.8 X IQ-^ujo- 
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FIG. 2: Emission spectra in triple dots for different initial 
conditions. The thin lines between the levels (in each dot) 
indicate the optical couplings between the dots with value 
G = 15.25. El ~ E2 = E2 ~ E3 = 0.02. F = 6.8 x lO^^c^o- 



The above neat theorem is very powerful in the "har- 
monic engineering". It leads to many important con- 
sequences and is very useful for the applications in de- 
signing the optical emission patterns of coupled nanos- 
tructures. In the following, we give a few examples of 
the applications of our theory. In the numerical calcula- 
tions, the density matrix is obtained by solving equation 
([2]) through Runge-Kutta method with the time step of 
0.01/wo, total steps of 1500000 and appropriate initial 
conditions. Emission spectra are obtained by numeri- 
cally calculating dipole through the evolution of density 
matrix elements pij (t) . We set huJo as the unit of energy. 

Double dots We first consider the system of double 
dots with two levels. One can define Q as the time shift 



6 : t ^ t + T/2 combined with fii : ci — > — ci (here and 
in the following Cj,j — 1, ...,N, refers to the annihilation 
operator for state \j) in the dot j). It is easy to see that 
under such transformation H is invariant, P = G{c^C2 + 
cjci) is odd; and the initial condition P22 = 1 is invariant 
under transformation fli. So there is no even component 
and odd peaks appear only, as seen in Fig. (la) and 
observed in many previous literatures. As our theory 
points out, the emission patterns depend on the initial 
condition. We can apply our theory with ■ ci < — ci 
for the initial condition pn = P22 = 0.5. In this case P 
is even. Thus there is no odd harmonic and only even 
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FIG. 3: Emission spectra in triple dots with same initial con- 
dition P22 ~ 1.0, but with different structures. The thin lines 
between the levels (in each dot) indicate the optical couplings 
between the dots with value G = 5.0. (a) and (c) are the 
cases with chain structure, (b) and (d) are of loop structures. 
In (a) and (b), Ei - E2 = 0.01, E2 - Es = 0.02; in (c) and 
(d) El- E2 = E2- E-i = 0.02. r = 6.8 X 10"*a;o. 
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FIG. 4: Emission spectra in quadruple dots with different 
structures/initial conditions. The thin lines between the levels 
(in each dot) indicate the optical couplings between the dots 
with value G = 6.3. (a) and (b) are the cases with loop 
structure, (c) and (d) are of additional cross couplings. Ei — 
E2 = Ea - Ei = 0.008, E2 - Ea = 0.02. F = 6.8 x 10"**tJo. 



harmonics appear as seen in Fig. (lb). (Note that these 
peaks may appear as very small split from exact even 
components.) It also agrees with our previous analytical 
perturbation results. It is a very striking fact that even 
the peak for ujq (corresponding to Rayleigh peak, which 
is often very pronounced in the usual emission spectra) 
disappears due to the particular symmetry. As observed 
in Fig. (Ic), there are both odd and even components for 
the initial condition pn = 0.3, P22 = 0.7. This pattern is 



related to the fact that both symmetries ■ 6 and • 9 
are broken. 

Triple dots We first show we can use different ini- 
tial conditions to realize quite different emission patterns. 
These results are displayed in Fig. 2. It is seen that the 
emission spectra may show odd only, even only, or both 
even and odd components. In particular, the emission 
may be completely quenched under appropriate condi- 
tion. These results are based on our theorem. The sym- 
metry property related to ili : C2 — — C2 leads to the 
disappearance of even harmonics in Fig. (2a). Symme- 
try Q = fl2 ■ 0, ^2 ■ ci i — > C3, results in the spectrum 
without odd harmonics shown in Fig. 2(b). Quite in- 
terestingly, both fii and ^2 are satisfied in Fig. (2c), 
which leads to the complete quenching of emission. This 
interesting phenomenon is a direct consequence of sym- 
metry and is independent of the frequency of incident 
light. Thus it is different from the quenching of emission 
due to coherent trapping. The breaking of both sym- 
metries fii and ^2 leads to emission with odd and even 
components as shown in Fig. (2d). 

Figure 3 shows the dependence of emission patterns 
on the structure of nanosystems. Figs. (3a) and (3c) 
show systems with chain structures and Figs. (3b) and 
(3d) show systems with loop structures. Systems (c) and 
(d) are of symmetric energy levels, which is absent in 
systems (a) and (b). The different emission patterns are 
the direct consequences of the symmetries generated by 
ill for (a), ^2 for (d), ili and ft2 for (c) and none of fii 
or ^2 for (b). 

Quadruple dots Our theorem is very powerful and 
can be applied to more general/complex nanostructures. 
Here we give one example of the tailoring the emission 
patters of quadruple dots. Careful design of structures of 
CQDs, or initial conditions leads to interesting emission 
spectra as shown in Fig. 4. The generators of the related 
symmetries are: f^i : C2 — —€2,04 — >■ — C4 for (4a); ^2 ■ 
C2 < — > C3, ci < — > C4 for (4b) and (4d); None for (4c). 

Our theorem on the essential roles of symmetry for the 
emission patterns emphasizes the dynamic symmetry of 
the system, instead of spatial symmetry. For a system 
with spatial inversion symmetry fJo: there may be odd 
harmonics only since the Hamiltonian is invariant under 
the symmetric operation Q ~ ■ 6 and P is odd under 
r^o- It is clear that dynamic symmetry is more basic. Our 
theory points out that the emission spectra depend not 
only on Floquet states (in particular their parity), but 
also on the initial state. Moreover we've found nanos- 
tructures with more general symmetries, which lead to 
more interesting emission patterns. Our studies suggest 
the effective methods for generation the emission spec- 
tra with odd harmonics only, even harmonics only, both 
odd and even harmonics, or even quenching of all compo- 
nents. While previous work studied only the generation 
of even harmonics in addition to the odd harmonics (i.e. 
both odd and even harmonics appear). Our theory can 
be applied to more complex configurations. 

We have considered the high-order harmonics genera- 
tion from symmetry point of view. In experiments, the 
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symmetric configuration could be realized by careful de- 
sign the nanostructures. For instance, the energy lev- 
els can be adjusted by changing the confining potential. 
They could also be tuned by applying appropriate gate 
voltages. The optical couplings among the dots could be 
tuned by changing the inter-dot barriers, or the polar- 
ization of the incident field. In some situation, one may 
need weaker conditions. For instance, one may obtain the 
emission spectra with even harmonics only by using the 
configuration of triple-dot (loop structure) with degener- 
ate levels for arbitrary inter-dot optical coupling. (It is 
a consequence of the symmetry generated hy Q = ■ 
with f2: Ci —Ci). Another issue is the preparation of 
the appropriate initial state. Some initial state, for ex- 
ample, that with electron staying in the lowest energy 
state, can be obtained easily. One may use some spe- 
cial designed laser pulse to prepare other types of initial 
states As also seen in section II, small derivation 

of the initial conditions may not change the picture. For 
instance, the main features of emission spectra remain 
unchanged if the derivation of occupation probability for 
the initial states is less than 5%. 



SUMMARY 

Based on the new theorem (on symmetry in the time 
domain and energy spectrum domain) we proved, we pro- 
vide methods to tailor photon emission patterns in driven 
nanostructures by tuning the symmetry of the system. 
Apart from the emission spectra with only odd harmon- 
ics, or both odd and even harmonics, we are able to ob- 
tain emission spectra with just even harmonics. In suit- 
able condition, the photon emission can even be fully 
quenched. Our methods for tailoring emission spectra 
apply for general coupled nanostructures. 
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